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Abstract. We present a new integrable partial differential equation found by 
Vladimir Novikov. Like the Camassa-Holm and Degasperis-Procesi equations, this 
new equation admits peaked soliton (peakon) solutions, but it has nonlinear terms 
that are cubic, rather than quadratic. We give a matrix Lax pair for V. Novikov's 
equation, and show how it is related by a reciprocal transformation to a negative flow 
in the Sawada-Kotera hierarchy. Infinitely many conserved quantities are found, as 
well as a bi-Hamiltonian structure. The latter is used to obtain the Hamiltonian form 
of the finite-dimensional system for the interaction of N peakons, and the two-body 
dynamics (N = 2) is explicitly integrated. Finally, all of this is compared with some 
analogous results for another cubic peakon derived by Zhijun Qiao. 
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1. Introduction 

The subject of this paper is the partial differential equation (PDE) 

u t ~ u xxt + 4:U 2 u x = 3uu x u xx + u 2 u xxx , (1) 

which was discovered very recently by Vladimir Novikov in a symmetry classification of 
nonlocal PDEs with cubic nonlinearity [18]. The perturbative symmetry approach [17] 
yields necessary conditions for a PDE to admit infinitely many symmetries. Using 
this approach, Novikov was able to isolate the equation (TjQ) and find its first few 
symmetries, and he subsequently found a scalar Lax pair for it, proving that the equation 
is integrable. Due to the u xxt term on the left hand side of ([1]), this equation is not an 
evolutionary PDE for u. However, taking the convolution with the Green's function 
g(x) = exp(— |x|)/2 for the Helmholtz operator (1 — d 2 ) gives the nonlocal (integro- 
differential) equation u t + u 2 u x + g * [3uu x u xx + 2(u x ) 3 + 3u 2 u x ] = 0. It is convenient to 
define a new dependent variable m to be the Helmholtz operator acting on u, in which 
case the equation (TjQ) can be more concisely written as 

m t + u 2 m x + 3uu x m = 0, m = u — u xx . (2) 

Henceforth we work with the above form of the equation. 

The work of Camassa and Holm [3], who derived the equation 

m t + um x + 2u x m = 0, m = u — u xx (3) 

from an asymptotic approximation to the Hamiltonian for the Green-Naghdi equations 
in shallow water theory, has attracted a lot of interest in the past fifteen years, for 
various reasons. To begin with, it is remarkable that the Camassa-Holm equation (J3]) 
approximates unidirectional fluid flow in Euler's equations at the next order beyond the 
KdV equation, and yet preserves the property of being integrable, fitting as it does into 
the hereditary symmetry framework of Fokas and Fuchssteiner [7J, with a bi-Hamiltonian 
structure and a Lax pair based on a linear spectral problem of second order. Also, 
while there are smooth soliton solutions of ([3]) on a non-zero constant background (or 
equivalently, with the addition of linear dispersion terms), the Camassa-Holm equation 
has peakon solutions, which are peaked solitons of the form 

N 

u (x,t) = ^2pj(t) exp(-|x - qj(t)\), (4) 
where the positions qj and amplitudes pj satisfy the system of ODEs 

N N 

Qj = J2 Pk e~ ]9j ~ 9hl , Pj = Pj ^Pk sgn(g J - q k ) e^^' 1 (5) 
fc=i fc=i 

forj = l,...,iV. 

The peakons are smooth solutions of ([3]) except at the peak positions x = qj, 
where the derivative of u is discontinuous. The equations ([5]) form an integrable finite- 
dimensional Hamiltonian system, corresponding to geodesic flow on an iV-dimensional 
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manifold with inverse metric gi k = exp(— \qj — qj,\). The positions qj and momenta 
Pj satisfy the canonical Poisson bracket {qj,Pk} = 8jk- The dynamics of two peakons 
(N = 2) was solved explicitly in the original paper by Camassa and Holm [3], while 
the explicit solution for arbitrary iV was found by Beals, Sattinger and Szmigielski [2]. 
Fuchssteiner also showed that the equation ([3]) is related via a reciprocal transformation 
to the first negative flow in the hierarchy of the Korteweg-de Vries equation. 

One might wonder whether the Camassa-Holm equation is the only integrable PDE 
of its kind, being a shallow water equation whose dispersionless version has weak soliton 
solutions. This turns out not to be the case. Degasperis and Procesi used an asymptotic 
integrability approach to isolate integrable third order equations, and discovered a new 
equation with the dispersionless form 

m t + um x + 3u x m = 0, m = u — u xx . (6) 

The Degasperis-Procesi equation turns out to be integrable, with a bi-Hamiltonian 
structure and a Lax pair based on a third order spectral problem [I], and it also arises 
in shallow water theory [6]. The equation ([6]) is related by a reciprocal transformation 
to a negative flow in the hierarchy of the Kaup-Kupershmidt equation [U [13] , and it 
also has peakon solutions of the form (j4|) whose dynamics is described by an integrable 
finite-dimensional Hamiltonian system with a non-canonical Poisson bracket (see [14J, 
and section 4 below). The explicit solution of the iV-peakon dynamics was derived by 
Lundmark and Szmigielski [16]. There are at least two distinct integrable analogues 
of the Camassa-Holm equation in 2+1 dimensions [T2l [15], while the Euler-Poincare 
equation on the diffeomorphism group (EPDiff) provides a geometrical generalisation 
of the Camassa-Holm equation in arbitrary dimension [11] . and admits weak solutions 
with support on lower-dimensional submanifolds. Rosenau also found various PDEs 
with nonlinear dispersion, which have solutions with compact support [21], some of 
which are relatives of the Camassa-Holm equation [19]. 

In what follows we present a bi-Hamiltonian structure for the integrable hierarchy 
of PDEs of which ([2]) is a member, present a matrix Lax pair corresponding to a zero 
curvature representation for this equation, and show how it is related via a reciprocal 
transformation to a negative flow in the Sawada-Kotera hierarchy. We also present a 
system of Hamiltonian ODEs for the dynamics of peakon solutions of ([2]), and explicitly 
integrate the equations for the interaction of two peakons. Finally, we compare our 
results with analogous properties of the integrable PDE 

m t + (m(u 2 - ul)j =0, m = u- u xx , (7) 

which was recently obtained by Zhijun Qiao [20]. Qiao's equation was the original 
starting point for our study, since it has cubic (rather than quadratic) nonlinear terms, 
and this is what led us to ask Vladimir Novikov to seek other integrable equations of 
this kind. 
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2. Lax pair and reciprocal transformation 



The equation (j2J) arises as a zero curvature equation F t — G x + [F, G] 
the compatibility condition for the linear system 



where 



/o 



1 



F*, 
G*, 

mX 





0, this being 



(8) 



1 

mX 




/ 



G 



^ - u 2 mX 

2_ 

3A2 
u_ 
A 



- u 2 mX 



\ 



V -V I ^ + uu x J 

We found the linear system (jSJ) directly by applying the prolongation algebra method of 
Wahlquist and Estabrook (see [22], and also [8]), but the details of this derivation will 
not be given here^]. In any case, once a Lax pair is given one can use it to derive most 
of the important properties of an integrable PDE. 

The first important observation we wish to make about Vladimir Novikov's equation 
is that it is connected to a negative flow in the Sawada-Kotera hierarchy via a reciprocal 
transformation. Upon rewriting the PDE (J2J) in the form 

{m 2 % + (m 2 /V), = 0, (10) 

it is immediately clear that m 2//3 is a conserved density Since each of the equations ([3]) 
and ([6]) has a conserved density of the form m 1 / 6 , for 6 = 2,3 respectively, and these 
densities yield reciprocal transformations to negative flows in more familiar hierarchies, 
this suggests that we should define the new independent variables X and T by 

dX = m 2/3 dx - m 2/3 u 2 dt, dT = dt. (11) 

The closure condition d 2 X = for the exact one-form dX in the reciprocal 
transformation (fTTj) is just the conservation law ([TO]) . Transforming the time evolution 
of m in (j2J), together with the definition m = u — u xx , leads to the equations 
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where V = m 2/3 and W = um * . The evolution equation for 1/V in the new 
independent variables X, T is the reciprocal transformation of the equation (j2j). 
However, in order to recognise ( |T2|) as a member of the Sawada-Kotera hierarchy we 
need to apply the reciprocal transformation to the Lax pair. (For details of the Sawada- 
Kotera hierarchy and its extensions we refer the reader to [TO].) 

By writing the column vector ^ in components as = (ipi, ip2, " l l ) z) T , we can 
eliminate ip\ and ip3 from *& x = F^ to get a single scalar equation for ip = ip2, namely 



ipxxx - 2m x m Vxx - (m x 



,m 1 — 2(m x ) 2 m 2 + l)ip x — fn 2 X 2 %jj. 



(13) 



X V. Novikov told us that he earlier found a scalar Lax pair for the PDE |2]) based on a third order 
spectral problem, by applying a reciprocal transformation to a symmetry of fifth order. Any scalar linear 
problem for ([2]) should be equivalent to the matrix system |8|). possibly after a gauge transformation. 
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When the reciprocal transformation (II ip is used to transform the x derivatives as 
d x = Vdx, the equation ffT3]) becomes 

Vxx (V x ) 2 1 
^xxx + Uij x = \ 2 ^, with l/=— ^ + L^_-_ (14) 

so that the second equation in (TT21 has the form Wxx + + 1 = for the same 
potential U. The third order operator d x + £/<9x in (1141) is the standard Lax operator 
for the Sawada-Kotera hierarchy, and by transforming the t derivatives in ^ t = G^Sf 
according to d t = dx — W 2 dx we find that the T evolution of ip is given by 

^ = J? ( W ^xx ~ W x ^x) - V»- (15) 

After gauging ip by a factor of e 2T ^ 3x2 ^ to remove the final term above, and then 
replacing A 2 by A and setting <fi = — 3W, we see that (fill) and ( |T5l) are respectively 
equivalent to equations (2.25) and (2.26) in [13], and the compatibility requirement 
ipTxxx = ipxxxr for this pair of scalar equations gives two conditions, namely that 
Wxx + UW is independent of X, and Ut + 3Wx = 0. The latter two conditions follow 
from (JT2J) provided that U is given in terms of V as in ( fl4|) . 



3. Conserved densities and bi-Hamiltonian structure 

The Lax pair (jHJ) can be used to find infinitely many conserved densities for (J2J). Upon 
setting p = (log^a; in ffl3|) it is clear that p satisfies the equation 

Pxx + "ipPx + P 3 - 2m x m~ 1 (p x + p 2 ) + (rn(m~ x ) xx -ljp = m 2 \ 2 . (16) 

The corresponding t evolution of ip implies that p t = F x for some flux F, and so by 
expanding p in powers of A one finds coefficients that are conserved densities. The 
asymptotic expansion for A — > oo has p 3 ~ m 2 A 2 , so p ~ m 2 / 3 A 2//3 , which extends to 
an infinite series p ~ m 2 ^ 3 \ 2 ^ 3 + Yl'jLi Pj^~ 2 ^ 3 - The densities pj are all determined 
recursively from ( fl6i) as local functions of m; for example pi = m~ 5 ^ 3 m xx — |m _8//3 m 2 + 
3m~ 2 / 3 . An expansion in positive powers of A for A — > can consistently begin with 
p ~ —muX 2 , but one must solve a second order differential equation to obtain each 
subsequent term, which leads to increasingly nonlocal expressions in m and u. Since we 
know that ([2]) is reciprocally related to a negative Sawada-Kotera flow, it is natural to 
regard the pj as densities for Hamiltonians that generate a positive hierarchy of flows, 
with the expansion around A = producing Hamiltonian densities for negative flows. 

Having found these conserved densities, we require a pair of Hamiltonian operators 
Bi, B2 which are compatible (in the sense that B\ + £>2, or any linear combination of 
them, is Hamiltonian) and can be used to generate the hierarchy of flows that commute 
with (j2J). From earlier studies on the Camassa-Holm and Degasperis-Procesi equations 
|13[ [TJ], we know that all nonlocal operators of the form 

B = m 1 - 1/b D x m 1/b Gm 1/b D x m 1 ~ 1/b , (17) 
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with G = (c\D x + C2-D^)~ 1 for constants b,ci,C2, are Hamiltonian, and have Casimir 
J m l l b dx. In fact, the case 6 = 2 gives the third Hamiltonian structure for the Camassa- 
Holm equation, and 6 = 3 gives the second Hamiltonian structure for the Degasperis- 
Procesi equation. Since J m 2 ^ 3 dx is a conserved quantity for ([2]), this suggests we should 
consider the operator (|T71) with b = 3/2, and indeed we find that the equation can be 
written in Hamiltonian form as 

m t = B\- — , H—-f mudx, (18) 
dm 4 J 

for £>i = — 18£>|& =3/ /2 in the case c\ = 4, c 2 = —1. Some other conserved quantities are 
H x = / \ (u 4 + 2u 2 u 2 x - f ) rfx, if 5 = Jm 2 / 3 cix, # 7 = J ±(m- 8 / 3 m 2 + 9m" 2 / 3 )rfx = 

/ fiidx, and the next one has leading term H n = f (Tn~~ m / 3 m xxx + . . .)dx (up to 
rescaling). These are the first few in the sequence of Hamiltonians that generate local 
symmetries of weight k = ±1 mod 6 according to 
„ 8H k SH k+6 

m ^ = ^ = Sl ^T' (19) 

where B 2 = (1 — D'fym^ 1 D x m~ l {\ — D 2 ). The recursion operator is 1Z = B^B^ , and 
it generates the flows lZ n m x of weight 6n + 1 and the flows lZ n m t5 of weight 6n + 5. 
However, when k = 5 or 7 the rightmost part of the identity f|T9l) fails, since both i7 5 
and i?7 are Casimirs for B±; and the Hamiltonian if is a Casimir for £>2. The proof of 
the following theorem will be presented in a forthcoming article. 

Theorem 1 The operators B\ = —2(3mD x + 2m x )(AD x — D 3 )~ 1 (3mD x + m x ) and 
Bi = (1 — D^mT 1 D x mT x {\ — D 2 ) provide a bi- Hamiltonian structure for the hierarchy 
of symmetries of the equation (0|). 

4. Peakon solutions 

From ffTUl) the travelling waves u = u(z), z = x — ct of ([2]) satisfy [u 2 — c)m 2 ^ = const. In 
the general case this gives m = ^c 2 D(u 2 — c)~ 3 / 2 for constant D / 0, which integrates 
further to (u') 2 = u 2 + cDu(u 2 — c) -1 / 2 + cE, for another constant E. This can be 
reduced to a quadrature which is the sum of elliptic integrals of the third kind, namely 

dz= , W w+lJ =, w = u(u 2 - c)- 1 / 2 . (20) 
2^{Dw + E){w 2 - l) + w 2 

However, if we require waves that vanish at spatial infinity, then D — 0, which implies 
that m = whenever u 2 ^ c. No smooth solution can satisfy the latter requirement, 
but this observation suggests that there should be a weak solution of the form 

u(x, t) = iVce" 1 ^" 1 " 1 , c> 0, a; constant, (21) 

which has the same form as the peakon for the Camassa-Holm and Degasperis-Procesi 
equations, except that the amplitude is the square root of the speed rather than being 
equal to the speed, as is the case for the peakon solutions of (j3]) and The expression 
( |2~TT) has m = away from the peak, and u 2 = c at the peak, but to regard it as a weak 
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solution of d2J) it is necessary to substitute it into the equation and integrate against 
suitable test functions with support around the peak. For the single peakon (121]) we 
have m = ±2y/c5(x — ct — x ), but there is some subtlety in interpreting this as a 
solution, because u x = ^^/csga(x — ct — xo)e~\ x ~ ct ~ x °\ and m are distributions, while 
the equation (j2J) includes the product u x m. The integrals can be regularised by taking 
the convention sgn(O) = 0, but a more rigorous alternative is to construct the peakon 
distribution as a limit of smooth solutions of the PDE. For the Camassa-Holm equation 
it is known that the single peakon arises as a weak solution in this way (see [TSj for 
a very detailed treatment), and multi-peakons arise similarly as a degenerate limit of 
algebro-geometric solutions pp. 

If we take u to be a linear superposition of N peakons, as in (j3J), so that 
m = 2^2j =1 Pj{t)S(x — qj(t)), then substituting into the equation (j2j) and integrating 
against test functions supported at x = qj gives the equations of motion for the peak 
positions and amplitudes. 

Proposition 1 The equation (Hj) has peakon solutions of the form Q), whose positions 
qj{t) and amplitudes pj(t) evolve according to the dynamical system 

Qj = Ejwftft^-* 1 - 1 *-* 1 , (22) 

The above equations are not in canonical Hamiltonian form. However, in [H] one of us 
showed how Hamiltonian operators of the form ( fl7|) are reduced to Poisson structures 
on the finite-dimensional submanifold of N peaks or pulses, resulting in the Poisson 
bracket 

{qj,qk} = G(qj - q k ), {qj,p k } = Q> - l)G\qj - q k )p k , ^ 
{Pj,Pk} = -(b-l) 2 G"(q j - qk)PjPk, _> 

where G is the skew-symmetric Green's function for the operator G. For N > 2, the 
Jacobi identity holds for this bracket if and only if G satisfies the functional equation 

G'(a)(G(/3) + G(^))+ cyclic = for a + /? + 7 = 0. (24) 

This functional equation is also a sufficient condition for the operator ( TXT?) to be 
Hamiltonian, and Braden and Byatt-Smith proved in the appendix to [H] that the 
unique continuously differentiable, odd solution of equation ( |24l) is G(x) = v4sgn(x)(l — 
e~ B \ x \) for arbitrary constants A, B. Up to rescaling x, this is the Green's function for 
the operator G = (D x — D^)^ 1 (or G = D' 1 in the degenerate case B — ► oo). In the case 
at hand, the operator B in Theorem [TJ has G = (4:D X — Z^) -1 , and the Hamiltonian H 
reduces to a conserved quantity h for the equations of motion (1221) . which is quadratic 
in the amplitudes pj. 

Theorem 2 The equations |HP for the motion of N peakons in the PDE (fj|) are an 
Hamiltonian vector field 

q j = { qj ,h}, Pj = {pj,h} 
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for the Hamiltonian h = fc=iP?Pfc ex P( — l<?i ~Qk\) > with the Poisson bracket specified 
by 

{Qj,Qk} = sgn(g i -g fe )(l-e- 2 l*- ( "=l), Uj,Pk} = e~ 2 ^- qk \p k , 

{Pj>Pk} = sgn(g j - qk)e~ 2lq i~ qkl Pjp k . 

We conjecture that the equations ( 1221) constitute a Liouville integrable Hamiltonian 
system with N degrees of freedom. For N — 1 this is trivial, and for N = 2 the result 
follows from the existence of a second independent integral in involution with h, namely 

fc=p^(l-e- 2191 - 921 ), {k,h} = 0. (26) 

The invariant k is degree four in the amplitudes, and for all N there is an analogous 
integral, quartic in pj, obtained by restricting the Hamiltonian Hi to the peakon 
submanifold. Indeed, the conserved densities for the negative flows in the hierarchy of 
the PDE (121) should all reduce to integrals for the iV-peakon dynamics, but the explicit 
construction of iV independent Poisson-commuting integrals for (1221) is still in progress. 
It is also worth mentioning that the Lax pair ([S]) can be used to obtain an N x N Lax 
matrix for the finite-dimensional system, satisfying 

L$ = -A 2 $, L = SPEP, (27) 

where S jk = sgn(^- - q k ), P = diag(pi, . . . ,p N ), E jk = exp(-|g j - q k \). The jth 
component of the vector <fr is just ip2(qj(t),t), where ip 2 (x,t) is the second component 
of \I/ in (jSJ), and the corresponding time evolution $ = M<3> yields the Lax equation 
L = [M, L] for the system (1221) . However, unfortunately the spectral invariants of L, 
which are the coefficients of the characteristic polynomial det(L + A~ 2 I) (a polynomial 
in A -2 ), do not provide enough integrals. For instance, when N = 2 we find that the 
trace of L vanishes, while the trace of L 2 gives k, but h does not appear. For higher 
values of N we have found that the spectral invariants of L have degrees 4, 8, 12, . . . 
but the integrals of degrees 2, 6, 10, . . . are missing. This leads us to expect that there 
should be another Lax representation for this system which would provide the correct 
number of integrals for Liouville's theorem. 

For the two-peakon dynamics, the equations of motion are 

qi = (pi + p 2 e-\ qi - q ^) 2 q 2 = (P2 + Pie- 1 " 21 ) 2 

p x = sgn(gi- q 2 )e- lqi - q2l (Pi+P2e- lqi - q2l )piP 2 , (28) 
p 2 = -sgn(<7i - q 2 )e~\ qi ~ q2 \p 2 + Pie~ lqi ~ q ^)p 1 p 2 , 

and without loss of generality we consider the case where the peaks are initially well 
separated, so that q\ « q 2 with q x ~ c%t, q 2 ~ c 2 t (for c\ > c 2 > 0), and we assume 
that both amplitudes are positive, so p\ — ► ^fc{ and p 2 — > tJo2~ as t — > — oo. In terms of 
these asymptotic speeds the Hamiltonian is h — ^(ci + c 2 ) and the quartic invariant is 
k = C\c 2 . Upon integrating the equations (1281) we find elementary formulae for p\ — p 2 , 
Pip 2 and e - ' 91-92 ', leading to the expressions 

P 2 -p\ = (ci - c 2 ) tanhT p x p 2 = J c x c 2 + ^g^sech 4 !' 

Q2 ~ qi = \ log (l + 1 % c 1 tff 2)2 cosh 4 T) , 
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where T = (c\ — c 2 )(t — t )/2, with t being an arbitrary constant. The formula for 
qi + q2 is somewhat more formidable, being given in terms of a certain quadrature as 

qi + Q2 = (ci + c 2 )(t - t ) + J f(T)dT + const. (30) 
The integrand / is 

2(c? - ci) f (ci - c 2 ) 2 + 8cic 2 cosh 2 T) 

/(T) = - j-^-, 

(ci - c 2 ) 4 + 16cic 2 (ci + c 2 ) 2 cosh 4 T 

and the quadrature can be performed explicitly by partial fractions in tanh(T), but the 
answer is omitted here. From (1291) and (1301) it is apparent that the peakons exchange 
speeds under the interaction, without a head-on collision, so that q\ ~ c 2 t, g 2 ~ c\t as 
£ — > oo. They also undergo a phase shift, which is described by the asymptotics of the 
term J f(T)dT in fl3"U|) . but the precise formula is rather unwieldy and will be presented 
elsewhere. 

5. Qiao's equation 

As we already mentioned, our interest in peakon equations with cubic nonlinearity began 
with Qiao's equation (jTJ), which can also be written as 

m t + (u — u x )m x + 2u x m = 0. (31) 

Qiao presented a 2 x 2 Lax pair for this equation given by the linear system ty x = U^, 
\Jr t = V* with 

U = 



V 




\(u 2 -u 2 x ) -A 1 (u - u x ) - \m\(u 2 - u 2 x ) 

u + u x ) + \m\{u 2 — u 2 x ) — A~ 2 — \(u 2 — u 2 ) 



2 """V"- ">xj " 2 

Qiao also found a bi-Hamiltonian structure for his equation, namely 

m t = B 1 — = B 2 - — (33) 
dm dm 

where 

B x = -W x mD x l mD x1 B 2 = -2(D X - D 3 X ), (34) 

and H, Hi are the same as the conserved quantities for ([2]) given in section 3 above. 
(In Qiao's original papers the quantity H , proportional to H here, is out by a factor 
of 2, while the quantity denoted Hi in [20j is missing the u% term.) Note that the 
first operator in (1341) is of the form (1171) with 6=1, and the compatibility of these 
Hamiltonian structures can be proved by a slight extension of a result in |13j . 
If we apply the reciprocal transformation 

Tft I 

dX = -^dx — -m{u 2 — u 2 x )dt, dT = dt 
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to Qiao's equation (J7|) then we find the pair of equations 

(m' 2 ) T = -2u x , (mux)x = ^{u/m - 1). (35) 

By transforming the Lax pair given by ( 1321) and writing a scalar linear problem for ipi, 
the first component of we find that the X part is 

ipi,xx + Ox - v 2 )^ = -\ 2 ip, v = m' 1 , 

which is the Schrodinger equation corresponding to the spectral problem for KdV, and 
the expression vx~ v 2 is the standard Miura map from modified KdV. The corresponding 
time evolution is 




a = u — mux/2, 



from which it is clear that the pair of equations fl35l) corresponds to a negative flow in 
the (modified) KdV hierarchy. 

Qiao has noted that the equation (J7J) does not have standard peakons of the form 
u = ce~\ x ~ ct \. The general travelling wave solution for this equation can be solved in 
terms of an elliptic integral, and some interesting wave shapes have been found in [20] in 
cases where this integral reduces to expressions in hyperbolic functions. However, here 
we should like to point out that, at least formally, peakons of the form u = ±y/ce~^ x ~ ct ^ 
(just as found for ([2]) above) do provide solutions of Qiao's equation. From the equation 
in the form (13TI) it is clear that if m is given by a delta function then the m 2 terms do not 
make sense. However, if we take travelling waves u = u(z), z = x — ct and integrate (171) 
along the z axis against an arbitrary test function <p, and then perform an integration 
by parts, we find 

J m(u 2 - {u'f - cj ip'(z) dz = 0. (36) 

For the peakon u(z) = y/ce~^ we have u'(z) = — y / csgn(2;)e _ ' 2: ' and m(z) = 2y / c<5(z), 
and this satisfies ( |36l) as long as we assume the usual convention that sgn(0) = 0. A 
more careful derivation could be carried out along the lines of [TH]. The equations for 
TV peakons should be extremely degenerate, since 6 = 1 and G(x) is proportional to 
sgn(a;) in the bracket ( 1231) . so pj are constant and the amplitudes of the peakons do not 
change. The same conclusion is reached by integrating (J7J against a test function and 
performing integration by parts. 

It seems that peakon equations with cubic nonlinearity have several novel features 
compared with the Camassa-Holm and Degasperis-Procesi equations, and there are 
many more things to be revealed by further study. 
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